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We show how a conditional displacement of the vibrational 
mode of trapped ions can be used to simulate nonlinear col- 
lective and interacting spin systems including nonlinear tops 
and Ising models ( a universal two qubit gate), independent 
of the vibrational state of the ion. Thus cooling to the vibra- 
tional ground state is unnecessary provided the heating rate 
is not too large. 
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One of the paths leading to the current interest in 
quantum computation begins with attempts to answer 
Feynman's question [0: can quantum physics be effi- 
ciently simulated on a classical computer? It is gener- 
ally believed that the answer is no, although there is no 
explicit proof of this conjecture. It then follows that a 
computer operating entirely by quantum means could be 
more efficient than a classical computer. Our belief in 
this conjecture stems from a number of algorithms, such 
as Shor's factorisation algorithm [||, which appear to be 
substantially (even exponentially) more efficient than the 
classical algorithms. A number of schemes have now been 
proposed for a quantum computer, and some have been 
implemented in a very limited way. What kinds of simu- 
lations might these schemes enable? A number of inves- 
tigators have attempted to answer this question . In 
this paper we consider this question in the context of the 
ion trap quantum computer model and show that there is 
a class of nonlinear collective and interacting spin models 
that can be simulated with current technology. 

Nonlinear collective and interacting spin models have 
long endured as tractable, nonlinear quantum models 
with wide ranging relevance. Such models have appeared 
in nuclear physics ||], laser physics [||, condensed mat- 
ter physics and of course as a theoretical laboratory 
to investigate aspects of nonlinear field theories . In 
many cases however the match between model and ex- 
periment is only qualitative. In this paper I will show 
how some of these models may be directly simulated on 
a linear ion trap with individual ion addressing as in the 
quantum computation architecture. 

The interaction Hamiltonian for N ions interacting 
with the centre of mass vibrational mode can be con- 
trolled by using different kinds of Raman laser pulses. 
A considerable variety of interactions has already been 



interaction that does not change the vibrational mode of 
the ions. Each ion is assumeds to be driven by a resonant 
laser field which couples two states, the ground state \g) 
and an excited state |e). The interaction Hamiltonian is 



Hi 



N 



(1) 



where Vli is the effective Rabi frequency at the i'th ion 
and we have assumed the dipole and rotating wave ap- 
proximation as usual. The raising and lowering operators 
for each ion are defined by ct- = \g){e\ and a+ = \e){g\. 
If we now assume that each ion is driven by an identical 
field and chose the phase appropriately, the interaction 
may be written as 



Hi = nnjy 



(2) 



where we have used the definition of the collective spin 
operators. 



where a — x,y, z and 
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achieved or proposed |12-14|. Consider first the simplest 



The interaction Hamiltonian in Eq |^ corresponds to a 
single collective spin of value j = N/2 precessing around 
the Jy direction due to an applied field. By choosing 
the driving field on each ion to be the same we have 
imposed a permutation symmetry in the ions reducing 
the dimension of the Hilbert space from 2^ to 2 + 1 . 
The eigenstates of may be taken as a basis in this 
reduced Hilbert space. In ion trap quantum computers 
it is more usual to designate the electronic states with a 
binary number as \g) ~ |0), |e) = The product basis 
for all N ions is then specified by a single binary string, or 
the corresponding integer code if the ions can be ordered. 
Each eigenstate, |j, to)^, of Jz is a degenerate eigenstate 
of the Hamming weight operator ( the sum of the num- 
ber of ones in a string) on the binary strings labelling the 
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product basis states in the 2 dimensional Hilbert space 
of all possible binary strings of length N. Collective spin 
models of this kind were considered many decades ago 
in quantum optics |^ and are sometimes called Dicke 
models after the early work on superradiance of Dicke 
| p5[ . In much of that work however the collective spin 
underwent an irreversible decay. In the case of an ion 
trap model however we can neglect such decays due to 
the long lifetimes of the excited states. However when 
the electronic and vibrational motion is coupled heating 
of the vibrational centre-of-mass mode can induce irre- 
versible dynamics in the collective spin variables. 

The natural variable to measure is Jz as a direct de- 
termination of the state of each ion via shelving tech- 
niques will give such a measurement. These measure- 
ments are highly efficient, approaching ideal projective 
measurements. The result of the measurement is a bi- 
nary string which is an eigenstate of J^. Repeating such 
measurements it is possible to construct the distribution 
for Jz and corresponding averages. Other components 
may also be measured by first using a collective rotation 
of the state of the ions. 

We now show how to realise nonlinear Hamiltonians 
using N trapped ions. By appropriate choice of Raman 
lasers it is possible to realise the conditional displacement 
operator for the i'th ion | |T^ , |l^ 

H = -ih{aia^ - a*a)ai''> (7) 

If the ion is in the excited (ground) state this Hamiltonian 
displaces the vibrational mode by a complex amplitude a 
(—a). In the case of N ions with each driven by identical 
Raman lasers, the total Hamiltonian is 

H — —ih(aa^ — a*a)Jz (8) 

By an appropriate choice of Raman laser pulse phases 
we can then implement the following sequence of unitary 
transformations 

Unl = e'''^"^''' eJ'^''^'^' e~^'^^^'^' e^'^''^'-'' (9) 

where X = {a + at)/V2, P = -i{a - a'')/V2. Noting 
that 

it is easy to sec that 

Unl = e-'"" (11) 

where = K^Kp which is the unitary transformation gen- 
erated by a nonlinear top Hamiltonian describing preces- 
sion around the Jz axis at a rate dependant on the z 
component of angular momentum. Such nonlinear tops 
have appeared in collective nuclear models and form 
the basis of a well known quantum chaotic system |ll9|] . 



It should be noted that the transformation in Eq([Tl|) 
contains no operators that act on the vibrational state. 
It is thus completely independent of the vibrational state 
and it does not matter if the vibrational state is cooled 
to the ground state or not. However Eq(|Tl|) only holds if 
the heating of the vibrational mode can be neglected over 
the time it takes to apply the conditional displacement 
operators. We discuss below what this implies for current 
experiments. 

In itself the unitary transformation in Eq ( pi] ) can gen- 
erate interesting states. For example if we begin with all 
the ions in the ground state so that the collective spin 
state is initially \j, —j)z and apply laser pulses to each 
electronic transition according to the Hamiltonian in Eq 
@) for a time T such that VLT — -k /2 the collective spin 
state is just the Jx eigenstate |j, —j)x- If we now apply 
the nonlinear unitary transformation in Eq ( |ll|) so that 
9 = 7r/2 we find that the system evolves to the highly 
entangled state 

1+) = ^(e--/^|j,~j% + (-l)-''e"/1j,j%) (12) 

Such states have been considered by Bollinger et al. 
Jl?! in the context of high precision frequency measure- 
ments, and also by Sanders | p^ . They exhibit interference 
fringes for measurements of Jz- As noted above a mea- 
surement of Jz is easily made simply by reading out the 
state of each ion using highly efficient fluorescence shelv- 
ing techniques. This particular nonlinear model is a well 
known system for studying quantum chaos, as we now 
discuss. 

The nonlinear top model was introduced by Haake 
|]l9| , ^ as a system that could exhibit chaos in the clas- 
sical limit on a compact phase space, and which could 
be treated quantum mechanically with a finite Hilbert 
space. This removed the necessity of truncating the 
Hilbert space and the possibility of thereby introducing 
spurious quantum features. The nonlinear top is defined 
by the collective spin Hamiltonian, 

H=^.Pz+pJy J2 S{t~nT), (13) 

where t is the duration between kicks, J = (Jx, Jy, Jz) 
is the angular momentum vector, and = j{j + 1) is a 
constant of the motion. 

As the Hamiltonian is time periodic the appropriate 
quantum description is via the Floquet operator 




which takes a state from just before one kick to just be- 
fore the next, i.e., \tp) — > U where Jz and Jy are the 
usual angular momentum operators, and j is the angu- 
lar momentum quantum number. The first exponential. 
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exp 



describes the precession about the 



2;-axis, and the second, Uk = cxp \^—i^Jyj , describes 
the kick. 

The classical dynamics can be reduced to a two dimen- 
sional map of points on a sphere of radius j , and the 
angular momentum vector can be parameterised in polar 
coordinates as 



J — j(sin 8 cos $, sin 8 sin $, cos 8). 



(15) 



The first term in the Hamiltonian ( |13[ ) describes a non- 
linear precession of the top about the z-axis, and the sec- 
ond term describes periodic kicks around the y-axis. The 
classical map for p — t:/2 and k = 3 has a mixed phase 
space with periodic elliptical fixed points and chaotic re- 
gions. 

It is now clear that this model can be simulated by 
the sequence of pulses in Eq ( ITTl) with appropriate values 
for the pulse area, together with a single linear rotation. 
This presents the possibility of directly testing a num- 
ber of ideas in the area of quantum chaos, particularly 
the idea of hypersensitivity to perturbation introduced 
by Schack and Caves |^ . Of particular interest here is 
the ability to very precisely simulate the measurement 
induced hypersensitivity discussed in |22j. In that pa- 
per the kicked top was subjected to a readout using a 
single spin that could be prepared in a variety of states. 
The interaction between the top and the readout spin is 
described by 



Ui — exp (^~i^Jy(j\, 



(R) 



(16) 



where we regard one ion as set aside to do the readout and 
label it with a superscript. It is relatively straight for- 
ward to generate this interaction via the pulse sequence 
of conditional phase shifts 



Unl = e 



iK.^Xai"'> iK„Pa["^ ^-iK^X.h jK„Pa["^ 



(17) 



with fi = K^Kp. It is now possible to consider a long 
sequence of measurements made at the end of each non- 
linear kick and record the resulting binary strings of mea- 
surement results. 

Initial states of the kicked top can be easily be prepared 
as coherent angular momentum states by appropriate lin- 
ear rotations. In the basis of orthonormal eigenstates, 
and \ j,m) = j{j + 1) \ j,m). the spin coherent states 
can be written as a rotation of the collective ground state 
p9|,p3[ through the spherical polar angles (0, (j)), 



I7) = exp i9{X sin </) - Jy cos </))) | j, -j) 



(18) 



where 7 = e'*tan {^) ■ This can be achieved by identi- 
cal, appropriately phased pulses on each ion separately. 
Initial states localised in either the regular or chaotic 



regions of the classical phase space may thus be easily 
prepared. 

Using a sequence of conditional displacement opera- 
tors that does distinguish different ions we can simulate 
various interacting spin models. As interacting spins are 
required for general quantum logic gates, these models 
may be seen as a way to perform quantum logical oper- 
ations without first cooling the ions to the ground state 
of some collective vibrational mode. 

Suppose for example we wish to simulate the interac- 
tion of two spins with the Hamiltonian 



Hint = ^xo-; 



(19) 



The required pulse sequence is 



Uint — e^'^"'^"''^^ e^'^^^'^'^^ e~^'^'''^"^^'' e^'^^'^"'^^ (20) 



This transformation may be used together with sin- 
gle spin rotations to simulate a two spin transformation 
that is one of the universal two qubit gates for quantum 
computation. For example the controlled phase shift op- 
eration 



JJ — ^j-i'^k)! (e|(»|e)2 {e| 
^ cp — ^ 



may be realised with x = tt as 



(21) 



(22) 



Once again this transformation does not depend on the 
vibrational state and so long as it is applied faster than 
the heating rate of the collective vibrational mode it can 
describe the effective interaction between two qubits in- 
dependent of the vibrational mode. 

We have proposed a scheme, based on conditional dis- 
placements of a collective vibrational mode, to simulate a 
variety of nonlinear spin models using a linear ion trap in 
the quantum computing architecture and which does not 
require that the collective vibrational mode be cooled to 
the ground state. However the scheme docs require that 
the heating of the collective vibrational mode is negligi- 
ble over the time of the application of the Raman condi- 
tional displacement pulses. It does not matter that the 
ion heats up between pulses. If the pulses were applied 
for times comparable to the heating times the pulse se- 
quences described above would not be defined by a prod- 
uct of unitary transformations but rather by the com- 
pletely positive maps which include the unitary part as 
well as the nonunitary heating part. Such maps provide 
a means to test various thermodynamic limits of nonlin- 
ear spin models and will be discussed in a future paper. 
In current experiments the heating time is estimated to 
be of the order of 1 ms, which is much shorter than the 
theoretically expected values that are as long as seconds 
]l2[. The source of this heating is unclear but efforts 
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are under way to eliminate it. so we can expect heat- 
ing times to eventually be sufficiently long to ignore. In 
current experiments however the sequence of conditional 
displacements would need to be applied on time scales of 
less than 1 ms. This is achievable using Raman pulses. 
We thus conclude that simple collective and interacting 
spin models with a few spins are within reach of current 
ion trap quantum computer experiments. 
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